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Abstract 

We study effects of charged Higgs boson exchange in the B semileptonic 
decays B -^ D^*>ti't-. Both branching ratio and r polarization are exam- 
ined. We use the recent experimental data on semileptonic B decays and the 
heavy quark effective theory in order to reduce theoretical uncertainty in the 



X 

;_( ■ hadronic form factors. Theoretical uncertainty in the branching ratio is found 



to be rather small and that in the r polarization is almost negligible. Their 
measurements will give nontrivial constraints on the charged Higgs sector. 
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The evidence for top quarks [|l| leaves the Higgs sector as the only missing part of 
the standard model (SM). In the minimal SM, we have one Higgs doublet, which gives a 
neutral scalar particle as a physical state. An extension of the Higgs sector is an interesting 
possibility for new physics beyond the SM. One of the most attractive possibilities is the 
supersymmetric extension of the SM [Q . In the minimal supersymmetric SM, we have to 
introduce two Higgs doublets in order to cancel the anomaly and to give the fermions masses. 
Another important possibility is CP violation in the Higgs sector [§,0. It is known that CP 
can be broken in the Higgs sector if we have two or more Higgs doublets. Since the existence 
of one or more charged Higgs bosons is an inevitable consequence of the multi-Higgs-doublet 
extensions of the SM, the search for their effects is one of the key points in the quest for 
new physics. 

The most stringent experimental bound on the charged Higgs boson mass at present is 
niH ~ 260GeV, given by the measurement of the inclusive radiative b quark decay 6 — *> 57 0. 
This bound was found for the two-Higgs-doublet model of the "SUSY-type" Higgs couplings 
to fermions, called Model II |^. Since this process takes place via 1-loop diagrams, the bound 
may be changed depending on details of the model considered. If the model contains new 
particles other than the charged Higgs boson which contribute to the b ^ s'-f process, the 
above lower bound may be modified. This is indeed the case in the minimal supersymmetric 
SM [^. From this point of view, it is worthwhile to investigate charged Higgs boson effects in 
tree level processes which are less dependent on the other sectors of the multi-Higgs-doublet 
models. 

In this paper, we study effects of the charged Higgs boson on the branching ratio and the 
r polarization of the processes B -^ D^*^tiJt-- These processes are expected to be much more 
sensitive to the charged Higgs sector than the semileptonic K decay processes because the 
Higgs couplings to fermions are proportional to the fermion mass. The 1% level branching 
ratio of these modes expected in the SM will give of order 10^ semi-tauonic B decay events 
at the planned B factories. 

In A^-Higgs-doublet models, we have A^ — 1 physical charged Higgs bosons. Their cou- 



plings to quarks and leptons are described by the following interaction Lagrangian [^: 

N-l 

Ch = {2V2GfY''' Y. [x^TJlVkmMdDr + YdJRMuVKMDL + Z.T^lMeEr] H^ + h.c. . 

(1) 

Here Hf is the i-th lightest physical charged Higgs boson, 

Ul(r) = {u, c, t)^(^) , Dl(r) = {d, s, b)l^R) , Nl = {ve, Jy,,, Vr)l , Er = (e, /i, r)^ , (2) 

Mu = diag.(m„, m^, m^) , Md = diag.(mrf, m^, m^) , Me = diag.(me, m^, m^) , (3) 

represent the quark and lepton fields and their masses respectively, and Vkm is the 
Kobayashi-Maskawa matrix. Note that the KM matrix which appears in the charged cur- 
rent mixing appears in the above Lagrangian. This is a consequence of the natural fiavor 
conservation ^ which we implicitly assumed in Eq. (|I]) in order to suppress fiavor changing 
neutral Higgs interactions. 

In the case of A^ > 3, the coefficients Xj, Yi, Zi can be complex, while they are real when 
N = 2. In particular they are real in Model II of two-Higgs-doublets or in the minimal 
supersymmetric SM, and can be written as 

Xi = Zi = tan/5 , Fi = cot /3 , (4) 

where tan/5 = Vu/vd is the ratio of the vacuum expectation values of the two Higgs fields. 

Given the above Lagrangian in Eq. (||) and the standard charged current Lagrangian, 
we can evaluate effects of the charged Higgs boson exchange in i? ^ MrVr processes for 
M = D OT D*. We adopt a helicity amplitude formalism since it is convenient for calculating 
r polarizations. We follow the convention in Refs. |]10| and p |. 

The W boson exchange amplitude is given by 

where \m = ±, denote three possible D* helicity states. Am = s stands for the D mode, 
and Ar = ± is the r helicity. The invariant mass squared of the leptonic system is g^, and 



X = Pr ■ Pb/^b is the r energy divided by the B meson mass in the B meson rest frame. 
The virtual W hehcity is denoted by Av^^ = ±, 0, s, and the metric factor rjxyy is given by 
rj±fl = 1 and rjs = (g^ — M^^)/M^r. The hadronic and leptonic amphtudes which describe 
the processes B -^ MW* and W* —>■ rvr are defined respectively by 



^A^(?') = C{Xw){M{pM. Am)|c7^(1 - i,)h\B{pB)) 



(6) 



and 



LtAq\x)=e^{\w){T{prAr)l^r{Pu)\fl^{l-l^)T^M 



(7) 



where e^{\w) is the polarization vector of the virtual W boson. Note that the leptonic 
amplitude L^^ depends on the frame in which the r helicity is defined. Ly^^ with the r 



helicity defined in the virtual W rest frame is given in Ref. [jTO| and the one with the r 



helicity defined in the initial B rest frame is given in Ref. [0. The hadronic amplitude is 
also given in Ref. [ p!0|] . 

The amplitude of the charged Higgs exchange can be written as 



MtM"^^)H = '^v^^Y.L'^ 



x,z: T"^' Hi^^' + Y,z: "^'^^^^ g^ 



« i\/r2 



Mfj^ - q 



2^^R 



* « i\/r2 



Mi - q 



2^^L 



where 



(8) 



r^M („2 



HZ{q') = 2(M(pM, \M)\cPR,Lh\B{pB)) , Pr 



1+75 



Pr 



1 -75 



(9) 



L^^{q\x) = {T{Pr,Xr)MP^)\r{^-l>r\0) . 



(10) 



Using the equations of motion, the hadronic and leptonic amplitudes of charged Higgs ex- 
change are related to those of W exchange with scalar W* polarization (A^y = s): 



H%L 



V¥ 



nib-mc 



Hs , H]^ =0 , if ^ 



VT rrO TT±fi _ TT±fl 



mil + irtc 



(11) 



■\t 






m^ 



(12) 



As can be seen in Eq. (pi]), the Higgs exchange does not contribute to the decay into 
transversely polarized D* meson {\m = ±) because of angular momentum conservation. 
Therefore, we study charged Higgs boson effects in B decays into D mesons {Xm = s) and 
those into longitudinally polarized D* mesons (Am = 0) in the B rest frame. 

Using the helicity amplitudes of Eqs. (^) and (^, it is straightforward to calculate the 
differential decay rateQ: 



^rA,, = 7^El-^tl'^'^3, (13) 



where -M^m ^ -^Aaj(?^' ^)w + M]^^^{q^ , x)h, and rf$3 = dq"^ dx/QAn^ is the three-body phase 
space. Also, r polarizations can be calculated conveniently with the helicity amplitudes. Let 
us consider the decay rate with a definite r spin direction. It can be written as 

1 r 



dTx.As) 



dTx,, + {dT';^,eL + dTi^e^ + dTl^^eT)-s , (14) 



where s is the unit vector which points toward the r spin direction in the r rest frame, and 
the basis vectors are defined as e^ = Pt/\Pt\, ^t = Vm x Pt/\Pm x Pt\, and e_L = e^ x e^, 
with the convention that the angle from pm to Pr lies between and tt. The situation is 
depicted in Fig. |l|. These definitions can be used both in the virtual W (or Higgs) rest frame 
and in the B rest frame. The components of the spin-dependent part of the decay rate in 
Eq. (|14D are given in terms of the helicity amplitudes as 

d^L = i^jMtJ-\M-,j)d^^^ (15) 

^rt, = ;^Re(^t.-^,-^.>^3, (16) 



^" TUB 



Uny{MtlMl^.)d^,. (17) 



^The decay distributions of the charge conjugate processes {B — > Mt^v^ with M = D, D*) are 
obtained by taking complex conjugate of all the couplings. 



From Eq. (0), the possible three r polarizations are defined by 

'' rfr(ei)+rfr(-ei) dV ' ^ ' 

_ dT{e^) - dT{-e^) _ dV^ 
"- dT{e±) + dT{-e±) dT ' ^ ^ 

_ dTjeT) - dTj-eT) dT^ 
^ ~ drier) + dT{~eT) dT ' ^^ 

where we omit the common helicity index Am- Note that Pi + P]_ + P^ = 1 at an arbi- 
trary kinematical configuration because of the definite neutrino helicity. The longitudinal 
polarization (Pl) and the perpendicular polarization (Pi) depend on the frame in which 
the r helicity is defined. On the other hand, the transverse polarization {Pt) is frame- 
independent. It is well-known that the transverse polarization is a T- violating quantity as 
long as the final state interaction can be ignored. The T- or CP-violating nature of the 
transverse polarization can be seen in Eq. (|17|) since all the tree-level amplitudes of Eqs. (^ 
and (|^) are chosen to be real in the CP-conserving limit in our convention. 

In order to find numerical predictions, the hadronic transition form factors are needed. In 
Ref . ||lO| , the following set of hadronic form factors are employed and the hadronic amplitude 
H)^'^' in Eq. (^) is given in terms of them: 

{D{pM)\cYb\B{pB)) = U{q'){pB+PMr + f4q'){PB - PmT , (21) 

{D*ipM,\M)\crb\BipB)) = ih{q^V"tlj,{pB+PM)p{pB-PM). , (22) 

{D*{pM, XM)\cYl,b\B{pB)) = f2{q')eZ + ^li ■ Pb {Mq'){pB + PmT + U{q'){PB - PmT} , 

(23) 

where eM = €m{pm, ^m) is the polarization vector of the D* meson. In the following, 
however, we adopt a different set of form factors which is more convenient to incorporate 



the results of the heavy quark effective theory. We use the following form factors |jT2[ : 



{D{v')\cYb\B{v)) = ,/^^^^^^[^4y){v + v'r + Uy){v-vr] , (24) 

{D*{v', XM)\n''b\B{v)) = i^J^^^^^iiv{yY''<"'t\,,v'^v„ , (25) 

(26) 

where v = Pb/i^b and v' = PM/f^Ai are the four-velocities of the B and M(= D, D*) mesons 
respectively, and y = v-v' = (m^ + m|/ — g^)/(2mBmM)- The form factors in Eqs. (^T])~(^3p 



can be written in terms of the form factors in Eqs. (^4|)^(^6D: 



/± = ±^[(l±r)e+-(lTr)^], (27) 

(28) 

where r = mM/rriB and we omitted the arguments of the /'s and ^'s. 

In the heavy quark limit and in the leading logarithmic approximation (LLA), we have 

ma 



e+ = ev = a = a3 = c7e, e- = a. = 0, (29) 

where ^(1) = 1 and C denotes the QCD correction factor in the LLA|^. We assume the 
following form of the universal form factor]^, 

((y^-ij^f, (30) 



^The factor C drops out in our results. 

^The following results in this paper are not affected significantly by adopting alternative forms 



such as those used in Ref. [|15| 



and determine the slope parameter p from the experimental data of semileptonic B decays 
||T6| . As a result, we obtain 

p= 1.08 ±0.11, (31) 

with X^irJ d.o.f. = 0.48. Eq. ( PT| ) gives the uncertainty of the predictions in the approxima- 
tion of Eqs. (IID and dO). 

In our numerical analysis, we consider only the effects of the lightest charged Higgs 
boson exchange. Other charged Higgs bosons (if they exist) are assumed to be too heavy to 
give significant contributions. Moreover, we concentrate on the case of the Model II Higgs 
couplings as shown in Eq. (|). So, the strength of the charged Higgs couplings to fermions 
are determined by tan/? only. 

We use mfe = 4.8GeV, rric = 1.4GeV, and m^ = 1.78GeV, and we do not consider the 
uncertainty in the quark masses, because its effect appears mostly through the combination 
mi,m^ ia,v? (3 / Mfj and a variation in mi, can be absorbed by changing tan/3 or the charged 
Higgs boson mass Mh. 

Our predictions for the branching ratio are shown in Figs. and ^. In Figs. |^(a) and ^(a), 
we show the decay rate of the process B -^ D^*^tiJt. in the presence of the charged Higgs 
boson exchange, normalized to the rate oi B ^ D^*^/^!?^ in the SM, against the charged 
Higgs mass for several values of tan (3. The shaded regions correspond to the predictions in 
the approximation of Eqs. (p9| ) and (|30|) within the uncertainty of Eq. (|3lD . In Figs. ||(b) and 
^(b), we also show the decay rate normalized to F, the decay rate of -B -^ D*^*^/iz/^ in the SM, 
but integrated in the same q^ region as the r mode, i.e. m'^ < i?^ < {tiib — niMY- As seen in 
Figs. and ^, this restriction in the q^ range decreases the uncertainty in the hadronic form 
factors and improve the sensitivity to the charged Higgs sector. The theoretical sensitivity 
to the charged Higgs sector can be represented by the minimum value oi R = M\y tan (3 /Mfj 
which can be detected in an ideal experiment. From Fig. 0(b), we expect the theoretical 
reach of i? ~ 6 with the uncertainty of Eq. (PT|). As can be seen in Fig. ^ the longitudinal 



D* {DD mode is less sensitive to the charged Higgs boson exchange because of the angular 
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momentum barrier. Actually, the hadronic amplitude H^{q'^) vanishes as q^ — > {niB — mMY- 

In order to get an idea on the size of the l/rrihc correction and the non-LLA QCD 

correction which violate the relation among the form factors in Eq. (|29|), we employ the 

estimation of these corrections by Neubert [|17|. In Ref. |jT^ both the l/nib^c correction as 



estimated by the QCD sum rule and the perturbative QCD correction beyond the LLA are 
given. The form factors in Eqs. (p^ ^ (|26|) are then written as 

Uy) = [«. + P^{y) + l^{y)] ^ , i = +,-,v,Ai, A^, A^ , (32) 

where a+ = ay = oai = c^As, = I, a- = Q.A2 = 0, /3i(|/) represents the perturbative 
QCD correction, and 7i(|/) is the l/m^.c correction. The functions /3j's and 7j's are given in 
Ref. [|l^. Assuming the form of ^{y) in Eq. (^) again, we obtain the following range of the 
slope parameter from the experimental data |T6|: 



p = 1.23 ± 0.09 , (33) 

with Xmin/d-o-f- = 0.55. 

The results on the decay rate by using Eq. (|32|) and the central value of Eq. (P3D, i.e. 
p = 1.23, are shown in Figs. ^ and |^ by dashed lines. The magnitude of the uncertainty 
from the range of the parameter p in Eq. ( ^3]) is roughly as the same as in the leading order 
approximation of Eqs. ( ^9]) ~ (|3T1) . Fig. ^ shows that the non- leading corrections are not the 



major uncertainty in the B to D mode, and the theoretical reach of -R ~ 6 for the D mode 
remains valid. On the other hand, as can be seen from Fig. |^, the non-leading corrections 
can be as large as the uncertainty in the slope parameter of Eq. (|33|) for the D^ mode. In 
future, the non-leading corrections may become dominant uncertainties in both the D and 
1^2 modes if the ranges of Eqs. (0) and ( |33D are reduced enough by detailed studies of the 
semileptonic B decays at B factories. 

As for the r polarizations, we can calculate any of their distributions by using Eqs. (|13|) 
and (|I^)~(^U|) given the hehcity amplitudes. For the couplings of Eq. (^, we obtain Pt = 0, 
and P| + Pf = 1 at any phase space point (g^,x). However, for simplicity, we concentrate 

9 



on the longitudinal polarization Pl integrated over the whole phase space separately in the 
numerator and the denominator of Eq. (|18|). Note that Pl + Pj 7^ 1 after this integration. 
In the following, we consider two of the possible frames in which the r helicity is defined, 
the virtual W{H) rest frame and the B rest frame. The longitudinal r polarization in the 
former frame is denoted by Pl{W*), and in the latter frame by Pl{B). 

Our numerical results on Pl{W*) and Pl{B) are given in Figs. § and |^. Fig. §(a) and 
(b) show the Pl{W*) and the Pl{B) respectively in the D mode. Fig. |^(a) and (b) show the 
same quantities in the D^ mode. The values in the leading order approximation of Eqs. (p9D 
and (|30|) with the uncertainty due to the range of p in Eq. (^l]) are again shown by the 



shaded region. The uncertainty in the prediction is found to be much smaller than that in 
the branching ratio; it is almost negligible except for Pl{B) in the D mode. The predictions 
with the non-leading corrections of Eq. (|32D for the central value of Eq. (|33D are also shown 
by dashed lines. The non-leading corrections can be regarded as the major uncertainty in 
the calculations of these polarizations except for Pl{B) in the D mode. 

From Fig. ^(a), we expect the best possible theoretical reach of -R ~ 4.5 among the 
calculations in this paper, regarding the difference between the center line of the shaded 
region which shrinks to almost a line and the dashed line as possible theoretical uncertainties. 
However, this estimation of the theoretical reach is rather ambiguous because it heavily relies 



on the estimation of the non-leading corrections in Ref. [T^. On the other hand, as expected, 
the D^ mode is less sensitive to the charged Higgs boson exchange in these polarizations 
too, see Fig. |. 

The above results on the sensitivity of the branching ratio and the r polarizations to 
the charged Higgs boson effects in the exclusive semi-tauonic B decays should be compared 
with those in the inclusive semi-tauonic B decay. The inclusive decay has been studied in 
Ref. W^, which finds the sensitivity of i? ~ 32 for the branching ratio, and i? ~ 20 for the 
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r polarization^, despite the fact that the uncertainties in the inclusive study of Ref. p8[ 



seem to be slightly smaller than those in the exclusive study of the present paper. These 
less sensitive results are understood as the above-explained insensitivity of the D* mode 
which gives a large portion of the inclusive decay and cannot be separated in the inclusive 
study. In other words, the D* modes dilute the effects of the charged Higgs boson exchange 
in the inclusive decay. On the other hand, in the exclusive study, we can select out the 
sensitive D mode. Therefore, measurements of the branching ratio and the r polarizations 
for B -^ DtUt mode may give the best bound on the charged Higgs boson exchange at the 
tree level. 

Finally, we comment on reconstruction of the r momentum. Reconstruction of the r 
momentum is desirable in several measurements discussed above, in particular it is neces- 
sary to measure x and g^ distributions. Even if we know the momenta pb, Pm, and ph in 
the decay process B -^ Mtv^ followed by r ^ hadrons(/i) + z/,-, the r momentum can- 
not be reconstructed because of two missing neutrinos. In this case, the r momentum is 
parametrized by its azimuthal angle in the virtual W{H) rest frame in which ph points to- 
ward the positive z direction. Improvements in vertex detector technology can improve the 
situation. Measurements of several quantities in semi-tauonic B decays should be improved 
significantly by systematically taking into account the vertex information. In principle, by 
using the knowledge about the tracks originating from M and r, we can measure the impact 
parameter between the flight lines of M and h. If the position of the B decay vertex or the r 
decay vertex is known in addition to this impact parameter, we can obtain the r momentum 
by determining the azimuthal angle mentioned above. However, there remains a two-fold 
ambiguity in general. To disentangle this ambiguity, both the B decay vertex and the r 
decay vertex should be known. In principle, the B decay vertex can be measured when the 
beam axis is well-known or when B decays into D* . The r decay vertex can be measured 



In Ref. |18|, the longitudinal r polarization in the B rest frame is discussed. 
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in the three-prong decays. 
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FIGURES 

FIG. 1. Kinematics and definition of the basis vectors of r polarization. 

FIG. 2. The branching ratios for the D mode: The shaded regions show the predictions within 
the uncertainty in p (Eq. (|3l|)) in the approximation of Eqs. ( P9|) and (pO|) , and the dashed hnes 
are the predictions with the non-leading corrections of Eq. (|3^ ) for the central value of Eq. (|33|). 
(a) The decay rate normalized to that of i^ ^ Dp.v^ in the SM: (b) The same as (a) except that 
the denominator is integrated over the region m^ < q^ < [piB — niMY- 

FIG. 3. The branching ratios for the D*^ mode: The same as Fig. ^ 

FIG. 4. The r polarizations for the D mode: The uncertainties are shown in the same way as 
Figs. |2| and y. The longitudinal r polarizations, (a) in the virtual W{H) rest frame, and (b) in the 
B rest frame are shown. 

FIG. 5. The r polarizations for the D*j^ mode: The same as Fig. ^. 
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